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Abstract 



Where are the degrees of freedom responsible for black hole entropy? 
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Considering the entanglement between quantum field degrees of freedom inside and outside the horizon 
as a plausible source of black hole entropy, we address the question: where are the degrees of freedom that 
give rise to this entropy located? When the field is in ground state, the black hole area law is obeyed and 
the degrees of freedom near the horizon contribute most to the entropy. However, for excited state, or a 
superposition of ground state and excited state, power-law corrections to the area law are obtained, and 
q-( more significant contributions from the degrees of freedom far from the horizon are shown 1 . 
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1 Introduction 

00 . 

0^ . The study of black holes (BHs) has always been a major testing arena for models of quantum gravity. 

The key issue has been to identify the microscopic origin of black hole entropy S BH . The questions that 
naturally arise in this context are the following: (i) Why is S BH , given by the well-known Bekenstein-Hawking 
relation [1,2], 
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S Bll = ~^ti {^p = V fcG/c 3 = Planck length, k B = Boltzmann constant^ (1) 

is proportional to the horizon area .A, as opposed to volume (usual for thermodynamic systems)? (ii) Are 
there corrections to this so-called 'area law' (AL) and if so how generic are these corrections? (hi) Can we 
locate the degrees of freedom (DoF) that are relevant for giving rise to the entropy? 

In the attempts to address these questions there have been two distinct approaches, viz., the one that 
associates 5 BH with fundamental DoF such as string, loop, etc. [3], and the other which attributes S Bli to 
the entanglement of quantum field DoF inside and outside the BH event horizon [4-7]. In this article we 
adopt the second approach and consider a quantum scalar field (in a pure state) propagating in the BH 
space-time. Since the BH horizon provides a boundary to an outside observer, the state restricted outside 
the horizon is mixed and leads to a non-zero entanglement (aka Von Neumann) entropy: S Ent = —k B Tr 
(plnp), where p is the mixed (or reduced) density matrix obtaining by tracing over the scalar DoF inside 
and outside the horizon. 

In Ref. [4, 5] — for scalar field is in the vacuum/ground state (GS) — it was shown that S Ent of scalar 
fields propagating in static BH and flat space-time (the DoF being traced inside a chosen closed surface) 
leads to the AL. In Refs. [8-10], the robustness of the AL is examined by considering non- vacuum states. 
It was shown that AL continues to hold for minimum uncertainty states like generic coherent state (GCS) 
or a class of squeezed states (SS), but for (1-particle) excited states (ES), or for a GS-ES superposition or 
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mixing (MS), one obtains power-law corrections to the AL. Although for large horizon area, the correction 
term is negligible, for small BHs the correction is significant. 

To understand physically the deviation from the AL for ES/MS, we ascertain the location of the micro- 
scopic DoF that lead to S Ent in these cases [10, 11]. We find that although the DoF close to the horizon 
contribute most to the total entropy, the contributions from the DOF that are far from the horizon are more 
significant for ES/MS than for the GS. Thus, the corrections to the AL may, in a way, be attributed to 
the far-away DoF. We also extend the flat space-time analysis done in [5] to (curved) spherically symmetric 
static black-hole space-times. 

In the next section, we first discuss the motivation for considering scalar fields for the entanglement 
entropy computations and then show that the scalar field Hamiltonian in a general BH space-time in Lemaitre 
coordinates, and at a fixed Lemaitre time, reduces to that in flat space-time. In sec. (3), we briefly review 
the procedure of obtaining the entanglement entropy and show the numerical results and estimations for 
the cases of GS, ES and MS. In sec. (4), we locate the scalar field degrees of freedom that are responsible 
for the entanglement entropy and compare the results for GS and ES/MS. We conclude with a summary 
and open questions in sec. (5). 

In the following, we use units with k B = c = H = 1 and set M% = 1/(16ttG). 



2 Hamiltonian of scalar fields in black-hole space-times 



Scalar fields can be motivated from the viewpoint of gravitational perturbations in static asymptotically 
flat spherically symmetric space-time background with metric g^ u . For a metric perturbation h^ u , the 
linearized form of the Einstein-Hilbert action is invariant under the infinitesimal gauge transformation 
hn V — > h^u + Imposing the harmonic gauge condition, i.e., d^{2h^ v — g^h") = [12] and keeping 

only the first derivatives of h^, one finally obtains the linearized spin-2 equation [13] 
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Assuming plane-wave propagation of the metric perturbations, i.e., h^ v = M pl e flu (p(x^) [e^ u = polarization 
tensor], in the weak-field limit, the above action reduces to the action for a massless scalar field <p propagating 
in the background metric g^ v : 
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Hence, by computing the entanglement entropy of the scalar fields, we obtain the entropy of the metric 
perturbations of the background space-time. The Hamiltonian of a scalar field propagating in a general 
spherically symmetric space-time background with line-element: 



ds 2 = -A( T ,£)dT 2 + 



is given by 
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where A,B,p are continuous, differentiable functions of (r,£) and we have decomposed ip in terms of the 
real spherical harmonics Z lm (6,4>), i.e., p{x^) = Y,i m L Pi m ( T ^) z im{0A)- 

In the time-dependent Lemaitre coordinates [13,14] the line-element is given by (4) with A(t,£) = 
1; B _1 (r, £) = 1 — f(r); p(r,£) = r(r,£). This line-element is related to that in the time-independent 
Schwarzschild coordinates by the following transformation relations [14]: 
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As opposed to the Schwarzschild coordinate, the Lemaitre coordinate is not singular at the horizon r n , and 
£(or, r) is space(or, time)-like everywhere while r(or, t) is space(or, time)-like only for r > r n . Choosing a 
fixed Lemaitre time (r = tq = 0), the relations (6) lead to: d^/dr = 1/^/1 — f(r). Plugging this in Eq.(5) 
and on performing the canonical transformations: II im — ► r\Jl — f(r) IL lm ; (p lm — > (p lm /r, the Hamiltonian 
reduces to that of a free scalar field propagating in flat space-time [15] 
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This holds for any fixed r, provided the scalar field is traced over either the region r G (0, r^] or the region 
r G [t7j,oo). Hence, evaluating the entanglement entropy of the scalar field in flat space-time corresponds 
to the evaluation of entropy of BH perturbations at a fixed Lemaitre time. 
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3 Entanglement entropy of scalar fields 

We discretize the scalar field Hamiltonian (7) on a radial lattice with spacing a: 
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where irimj are the momenta conjugate of <pi m j, (N + l)a is the infrared cut-off. H[ m in Eq.(8) is of the 
form of the Hamiltonian of N— coupled harmonic oscillators (HOs): 
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where the interaction matrix K^j is given by: 
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(10) 



The last two terms denote nearest-neighbour interactions and originate from the derivative term in (7). 
The most general eigenstate of the Hamiltonian (9) is a product of N— HO wave functions: 
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where x = Ux, (U T U = In), Kd = UKU T (diagonal), and v \ (i = 1 . . . N) are the indices of the Hermite 
polynomials (7i v ). The frequencies are ordered such that koi > kpj for i > j. 

The reduced density matrix is obtained by tracing over first n of the N oscillators: 



/n 
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It is not possible to obtain a closed form expression for p(x;x') for an arbitrary state (11). We resort to the 
following cases to compute the entropy numerically 2 using the relation S = Tr (plnp): 



2 The computations are done with a precision of 0.01% and for N = 300 and n = 100 — 200. 



(i) Ground state (GS) with A-particle wave-function: ip (x;x') ~ exp 

(ii) Excited (1-particle) state (ES) with A-particle wave-function: ipi(x; x') = V2a T Kp 2 x. ipo(x; x'), where 

a T = (ai, . . . , ajv) are the expansion coefficients and normalization of ip\ requires a T a = 1). We choose 
a T = l/y / o(0, • • • 0, 1 • • • 1) with the last o columns being non-zero. 

(hi) GS-ES linearly superposed (i.e. mixed) state (MS) with A-particle wave-function: ip(x; x') = 
CQipo(x;x') + c\ij)\{x;x'). Normalization of tp requires constants Co and c\ related by c\ + c\ = 1. 
For simplicity, we choose Co = 1/2. (See details in Ref. [10]). 

(a) In S vs In (R/a) (b) Asymptotic MS/ES entropy (c) Relative MS entropy 
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Figure 1: (a) log(entropy) vs log(i?/a) for GS, MS and ES (Eq/Hi). R = a{n + 1/2) is the horizon radius 
and N = 300, n = 100 - 200, o = 50. (b) Plots of S MS /S GS , S ES /S GS and S GS /S M s, S GS /S E s (for o = 50) 
with A to show the asymptotic nature of MS and ES entropies, (c) Best fit plots of Sms/S G s vs A for 
o = 30,40,50. 
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Results : For GS, one recovers the AL — S G s ~ A/a 2 , where a is the ultraviolet cutoff at the horizon (set 
to be ~ £p). For MS and ES, we obtain power-law corrections to the AL: 

Sms/es = S G s + cr{A/a 2 ) 1 ' y , (13) 

where a = constant of order unity and v is a fractional index which depends on the excitation o. As the 
horizon area A increases, the correction term becomes negligible and Sms &GS asymptotically. For 
small BHs, however, the correction is significant. Fig. 1 shows the logarithm of entropy versus log(i?/a) 
characteristics for GS, MS and ES (i? being the horizon radius), as well as the asymptotic equivalence of 
GS and MS/ES entropies, and the numerical fit that leads to the above result (13). 



4 Location of the degrees of freedom 

Let us take a closer look at the interaction matrix K{j, Eq.(10), for the system of HOs. The last two 
terms which signify the nearest-neighbour (NN) interaction between the oscillators, are solely responsible 
for the entanglement entropy of black holes. Let us perform the following operations: Operation I : We 
set NN interactions to zero (by hand) everywhere except in a 'window', such that the indices i,j run from 
q — s to q + s, where s < q. We thus restrict the thickness of the interaction region to t = 2s + 1 radial 
lattice points, while allowing it to move rigidly across from the origin to a point outside the horizon. The 
variation of the percentage contribution of the total entropy Stat for a fixed window size of t = 5 lattice 
points, i.e., pc(q) = [S(q,t = 5) /Stat] x 100, as a function of q is shown in Fig. 2 for N = 300, n = 100, in 
each of the cases GS and MS, ES with o = 30, 50. In all the cases pc(q) = when q is far away from n (i.e., 
horizon), whereas for values of q very close to n there are significant contributions to Stot- For GS, pc(q) 
peaks exactly at q = n. For MS and ES, however, the peaks shift towards a value q > n, and the amplitudes 
of the peaks also decrease as the amounts of excitation o increase [10,11]. Therefore (a) the near-horizon 
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Figure 2: Plots of the percentage contribution pc(q) to the total entropy as a function of window position 
q, for a window size t = 5 and fixed N = 300, n = 100, in each of cases of GS, MS and ES. For MS and ES 
the solid curve is for o = 30 whereas the broken curve is for o = 50. 
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Figure 3: The left panels show the variations of the percentage contribution pc(t) of total entropy with 
window width t for GS (o = 0) and ES (with o = 30, 50). The right panels show the plots of Apc(t) vs n — t 
for GS and ES. Both sets of panels are for N = 300 and n = 100, 150, 200. 



DoF contributes most to Stat, an d (b) the contributions from the far away DoF are more for MS and ES, 
than for GS. 

Operation II : We set the NN interactions to zero (by hand) everywhere except in a window whose 
center is fixed at p < i, j < n, and the window thickness t = n — p is varied from to n, i.e., from the 
origin to the horizon. For GS we find that about 85% of the total entropy is obtained within a width of 
just one lattice spacing, and within a width of t = 3 the entire GS entropy is recovered. Thus most of the 
GS entropy comes from the DoF very close to the horizon and a small part has its origin deeper inside. For 
ES, however, the total entropy is recovered for much higher values of t (than for GS) since the DoF that 
are away from the horizon contribute more and more as the excitation o increases. Thus larger deviation 
from the area law may be attributed to larger contribution to the total entropy from the DoF far from the 
horizon. The left panels of Fig. 3 depict the variation of the percentage contribution to the total entropy, 
i.e., pc(t) = [S(t)/Stot] x 100, as a function of t for GS (o = 0) and ES (with o = 30,50). The situation is 
intermediate for MS (which interpolates between the GS and the ES), i.e., the total entropy is recovered for 
values of t greater than that for GS, but less than that for ES (with same value of o). The percentage increase 
in entropy when the interaction region is incremented by one radial lattice point, Apc(t) = pc(t) —pc{t — 1), 
versus (n — t) plots for GS and ES are shown in the right panels of Fig. 3. In the case of GS the inclusion 
of the first lattice point just inside the horizon leads to an increase from to 85% of the total GS entropy. 
The next immediate points add more to this, but the contributions are lesser and lesser with inclusion of 
points further and further from the horizon. For ES however, inclusion of one lattice point adds 70(50)%, 
for o = 30(50), to the entropy, while the next immediate points contribute more than those for the GS. 



5 Conclusions 



We have thus shown that if the black hole entropy is looked upon as that due to the entanglement between 
scalar field DoF inside and outside the horizon, there are power-law corrections to the Bekenstein-Hawking 
area law when the field is in excited state or in a superposition of ground state and excited state. Although 
such corrections are negligible for semi-classical BHs, they become increasingly significant with decrease 
in horizon area as well as for increasing excitations. The deviation from the AL for ES and MS may be 
attributed to the fact that the scalar field DoF that are farther from the horizon contribute more to the 
total entropy in the ES/MS cases than in the case of GS. The near horizon DoF contribute most in any 
case, however. We have also extended the flat space-time analysis done in [5] to static spherically symmetric 
black hole space-times with non-degenerate horizons. 

We conclude with some open questions related to our work: (i) Can a temperature emerge in the 
entanglement entropy scenario and would it be consistent with the first law of BH thermodynamics? (ii) Is 
dS/dt > 0?, i.e., the second law of thermodynamics valid? (iii) Will the entanglement of scalar fields help 
us to understand the information loss problem? We hope to report on these in future. 
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